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"t^": ABSTRACT 
O 

o 

We explore the correlation between morphological types and mid-infrared 
(MIR) properties of an optically flux-limited sample of 154 galaxies from the 



Forth Data Release (DR4) of Sloan Digital Sky Survey (SDSS), cross-correlated 
with Spitzer SWIRE (Spitzer Wide- Area InfraRed Extragalactic Survey) fields 
of ELAIS-N1, ELAIS-N2 and Lockman Hole. Aperture photometry is performed 
on the SDSS and Spitzer images to obtain optical and MIR properties. The 
morphological classifications are given based on both visual inspection and bulge- 
disk decomposition on SDSS g— and r— band images. The average bulge-to- 
total ratio (B/T) is a smooth function over different morphological types. Both 
the 8yum(dust) and 24/im(dust) luminosities and their relative luminosity ratios 
to 3.6/im (MIR dust-to-star ratios) present obvious correlations with both the 
^ ■ Hubble T type and B/T. The early-type galaxies notably differ from the late- 

types in the MIR properties, especially in the MIR dust-to-star ratios. It is 
■ suggested that the MIR dust-to-star ratio of either vL u [8fj,m(dust)]/ vL v [3.6p,m] 

or vL„[24fim(dust)]/vL u [3.6iJ,m] is an effective way to separate the early-type 
galaxies from the late-type ones. Based on the tight correlation between the 
stellar mass and the 3.6/im luminosity, we have derived a formula to calculate 
the stellar mass from the latter. We have also investigated the MIR properties 
^ ■ of both edge-on galaxies and barred galaxies in our sample. Since they present 

similar MIR properties to the other sample galaxies, they do not influence the 
MIR properties obtained for the entire sample. 

Subject headings: galaxies: morphology — galaxies: formation — infrared: galax- 
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INTRODUCTION 



Ever since Hubble's famous paper outlined his classification system (jHubbldll926l . Il936l ) . 
morphological classification in conjunction with physical measurement has become an im- 
portant tool in extragalactic astrophysics. A number of quantitative classifiers have been 
developed or extended over the years to probe the struct ure of galaxies. There are para- 
metric classifi ers like radial mu l ti-Ga ussian deconvolution (jBendineUil Il99ll ) . bulge-disk de- 
composi tion (Bvun fc Freeman Il995 ). etc., and the nonparametric ones such as the C-A 
syste m ([Abraham et al .11 19941 . Il996l ). artificial neural nets train ed from visual classification 
sets ( jOdewahn et al.lll996l ). Gini Coefficient (ILotz et al.l 120041 ) and so on. As one of the 
main quantitative criteria and a function of the Hubble classification, the bulge-disk de- 
composition has now been widely u sed as an effective method to examine galaxy structures 
and morphologic al properties, (e.g.. 
Allen et al.lbood . etc). 



de Jondll996l : iBaggett et al.lll998l : iTasca fc Whitdl2005 



The variations of galaxy physical properties with m orphology and environment are cru- 



cial i n our understanding of the evolution of galaxies (IKennicuttl Il998l ; iBrinchmann et al. 



20041 ) . Numbers of properties such as the inte grated birthrate va r iation (ISandage 



the optical and infrared photometric pr operties (IBoselli et al.ll200ll ; IShimasaku et al. 



986) 



2001 



Popescu et al.ll2002l ; iDavoodi et al.ll2006l ). the star formation properties of galaxies in clu sters 
( jYuan et al. I l2005h . the circumnuclear Ha luminosity and bar structures (IShi et al.ll2006l ). etc, 
which have been investigated present regularity or correlation along different morphological 
types. Through substantial former investigations, it is known that as a result of different 
stellar populations and the amount of dust and gas for the environment of star-forming, the 
early-type galaxies (ellipticals and lenticulars) exhibit rather different properties compared 
with the late- type ones (spirals, irregulars and so on). 

As we know, averag ely one third of the total luminosity fr om normal galaxies is absorbed 
and re-radiated by dust (jMathislll990l ; |Popescu fc Tuffsll2002l ). and even higher fraction from 
galaxies with the most active star-forming activity (e.g. luminous infrared galaxies (LIGs), 
Sanders fc Mirabel! 1 19961 ; IWang et al.ll2006l ). B oth the Infrared Sp ace Observatory ( ISO - 



Kessler et al. 



19961 ) and Spitzer Space Telescope ( Werner et al. 2004 ) continued to explore the 



impo rtance of dust. Studies of dusty starburst galaxies (IPoggianti fc Wull2000l ; iFlores et al. 



20041 ) have shown that, most of the activities (e.g., star-formation and/or AGN emission) 
in these galaxies are hidden by dust, and the bolometric luminosities of the active systems 
are mostly emitted in the infrared. This suggests that the infrared emission is a sensitive 
tracer of the young stellar population and star formation rates (SFRs), and suffers weaker 
extinction. The MIR dust emissions mainly consist of polycyclic aromatic hydrocarbons 
(PAHs) emission and the continuum emission feature of warm dust component. Both of 
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these two c omponents have been investigated as reliable measures of the SFRs of galaxies 
as a whole (jPeeters et al.ll2004l ; IWu et al.l l2005al ) . Whereby studies in the MIR properties 
will certainly improve our understanding of the galaxies which have been well studied in the 
optical bands, and give an insight into the details of their star formation histories. 

This work tries to explore the relationships between the morphology and the MIR prop- 
erties, for a flu x-limited sample of normal gala xies which were selected from the galaxies 
of SDSS-DR4 ( Adelman-McCarthy et al.l 120061 ) cross-correlated with the Spitzer SWIRE 
( jLonsdale et al.ll2003l ). Considering that the Spitzer IRAC (Fazio et al. 2004) 8/xm band 
covers the strongest PAH feature (7.7/im), and the MIPS (Rieke et al. 2004) 24fim band 
covers the continuum emission of very small grains (VSGs) free of PAH features, we tried to 
employ the IRAC bands and MIPS 24fim band to investigate the MIR dust properties of en- 
tire galaxies. We performed elliptical aperture photometric analysis in both the optical and 
MIR bands. The 8/xm dust and 24/im dust luminosity and their dust-to-star ratios are used 
to quantitatively investigate the MIR properties of our sample galaxies. Two ways to clas- 
sify the morphologies of the sample galaxies were adopted: the T type of the revised Hubble 
classification system ( jde Vaucouleurs et al.lll976l ) by visual inspection, and the quantitative 
parameter of B/T by the bulge-disk decomposition. 

In §2, we describe the infrared and optical data, the sample construction, and the ellip- 
tical aperture photometry. The morphological classifications are presented in §3. In §4, we 
analyze the optical and MIR colors, and give the statistical results of the MIR dust properties 
against galaxy morphological types and the corresponding discussions. The conclusions are 
presented in §5. Throughout this paper, we assume a Hubble constant H = 70kms _1 Mpc _1 
and Qm = 0.3, = 0.7 in calculating the distance and the luminosity. 



DATA AND DATA REDUCTION 



2.1. The Sample 



We used the IRAC 3.6, 4.5, 5.8, 8.0fim and the MIPS 24/zm images from the northern 
SWIRE fields of Lockman Hole, ELAIS-N1, and ELAIS-N2. The BCD (Basic Calibrated 
Data) images of the four IRAC bands were obtained from Spitzer Science Ce nter, which 



inclu de flat-field corrections, dark subtraction, linearity and flux calibrations (jFazio et al. 



20041 ). The IRAC images (in all four bands) were mosaiced from the BCD images after 
pointing refinement, di stortion correction and c osmic-ray remova l with the final pixel scale 
of 0.6" as described bv iHuang etaD koo4 \ and lWu et aD J2005ah : likewise the MIPS 24/xm 
images were mosaiced in the similar way with the final pixel scale of 1.225" (jCao fc Wu 
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20071 : IWen et all 120071 ). Sources detected by SExtractor flBertin fc Arnoutslll996h in IRAC 
four bands and MIPS 24/im were matched with the Two Micron All Sky Survey (2MASS) 
sources to achieve astrometric uncertainties of around 0.1". 

The SDSS data provide full coverage of the SWIRE fields of Lockman and ELAIS-N2 
but cover only one third of ELAIS-N1 field. The ugriz corrected frames of spectroscopically 
observed galaxies were taken from the SDSS-DR4. The pixel scale of SDSS images is 0.4" 
( jStoughton et al.l 120021 ). The SDSS-DR4 ugriz spectrophotomet ric catalogue was cross- 
correlated with SWIRE MIR catalogue measured by SExtractor ( jBertin &: Arnoutsl Il996l ) 
by a radius of 2". The total survey area of the three northern SWIRE fields is ~24 deg 2 , of 
which the overlap wit h SDSS is ~15 deg 2 . To obtain the reliable morphological classification 
( IFukugita et al.ll2004l ). 163 bright galaxies with Petrosian magnitude r ^ 15.9 were selected. 
Only four SDSS galaxies were not matched with SWIRE MIR sources by 2", since the 
SDSS fiber observations mistakenly pointed to the off-nucleus regions rather than the nuclear 
regions of galaxies. This provides a preliminary sample of 159 galaxies. 

We excluded further five galaxies in our magnitude limited sample of 159 galaxies, 
because they failed in either the sky-background subtraction or the bulge-disk decomposition, 
or were severely contaminated by nearby bright stars. This led to a SDSS r— band flux- 
limited final sample of 154 galaxies. 142 of this sample have been imaged by IRAC four 
bands and 137 have been imaged by MIPS 24/im band. Finally, a sub-sample of 125 galaxies 
has images in all five MIR bands, and is used for further statistical discussion in §4.41 All of 
theses objects are local, with redshift less than 0.13. In our sample, six galaxies were with th e 
absolute B magnitude fainter than -18, and thus classified as dwarf galaxies (lMateolll9 98'l. 
Here , the B magnitude can be obtained from the SDSS g— and r— magnitudes (ISmith et al. 



20021 ). The distributions of SDSS r— band Petrosian magnitudes, the redshift, and B-band 
absolute magnitudes for the 154 sample galaxies as well as those for the 125 sub-sample 
galaxies are plotted in Figure [TJ All these distributions show that the sub-sample can well 
represent the flux-limited sample. 

Among the 154 sample galaxies, 133 with emission line detectio ns could be classified 



with the traditional log(Q I Il/H/3)-log(N II/Ha) diagnostic diagram ([Baldwin et al.l Il981 



Veilleux fc Osterbrockl Il987l ). The emission line fluxes were der ived from the catalogue of 
Tremonti et al.l (120041 ) . The criteria given by lKewley et al.l (120011 ) was adopted to distinguish 
the potential star-forming galaxies from AGNs, as is shown in Figure [2j 
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2.2. Photometry 



To obtain the accurate photometry, sky background fitting and subtraction were done 



(iZheng et al.lll999l ; lWu et al.ll2002l . l2005bl ) on both the SDSS corrected frames and the Spitzer 
images. All objects detected by SExtractor were masked to generate a background-only 
image, and the fitting sky-background was then subtracted . Photometry was performed on 
the background-subtracted images by IRAF task ELLIPSE (I Jedrzejewskilll987l ) . To embrace 
almost all the flux of these extended sources in the different bands, an elliptical isophote with 
the B-band surface brightness of 26 mag arcsec" 2 was adopted as the photometric aperture, 
based on the SDSS g— band images. With such elliptical apertures, the total fluxes were 
measured in all the wavelengths including the SDSS ugriz, the IRAC four bands and the 
MIPS 24/im band (see sample apertures marked in Figure [3]) with IRAF task POLYPHOT. 
Note that the point source functions (PSFs) of MIPS 24/zm are rather extended, therefore 
further aperture corrections were performed on this band. For objects with the equivalent 
radius of the elliptical apertures smaller than 15", aperture corrections were applied to 
calibrate the integrated flux to an equivalent radius of 15". The photometric accuracies 
of these different bands are quite s mall, less than 0.0 3 mag on average. Flux calibration 



accur acies of the IRAC four bands (IFazio et al.l 120041 ) and MIPS 24/zm band (IRieke et al. 
20041 ) are less than 10%. The final errors include both the above errors. 



The Galactic extinction in each SDSS filter from SDSS-DR4 cat alogue was ad opted, 

). The 
Kcor- 



and then intrinsic extinction was derived from the Balmer decrement (Calzctti 



photometric fT-correction was calculated using the method of iBlanton et al.l (12003 



200 



rect V4-1-4). No extinction correction has been performed on photometric results in the 
MIR bands since extinction effect is rather negligible in the infrared compared with optical 
wavelength. Due to the fact that our sample are all low redshift galaxies and, as of this 
work, there are no reliable .fT-correction for these MIR spectral ranges available yet, we ap- 
plied no i^-c orrection to the MI R photometry. All the measurements were converted to AB 
magnitudes (jOke fc Gunnlll983l ). 



Although PAH and VSG emissions dominate 8/im and 24/im bands for the majority 
of our sample, there is still a stellar continuum in these bands, especially for the early-type 
galaxies. Herewith a subtraction of the stellar contribution using the 3.6/xm luminosity 
was adopted , with a scale factor of 0.232 for the 8/xm band and 0.032 for the 24fim band 
(IHelou et al.l 120041 ) . Hereafter, we denoted the 8/xm(dust ) and 24/xm(dust ) representing the 
dust emissions after subtracting the stellar contribution (jWu et al.ll2005af ). 
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MORPHOLOGICAL CLASSIFICATION 



We conducted the morphological classification by tw o methods: visual inspection an d 
bulge-disk decomposition with GIM2D (Galaxy IMage 2D: ISimardlll998l ; ISimard et al.ll2002l ). 



3.1. GIM2D Fitting 



GIM2D is a two-dimensional photometric decomposition fitting algorithm which fits 
each image to a superposition of an bu lge componen t with a Sersic pro file, and a disk 
component with an exponential profile (jSimardl Il998l ; ISimard et al.l 120021 ) . GIM2D was 



employed to obtain the structural param eters of gala xies in our sample. The bulge component 
of the model is a profile of Sersic form (ISersid Il968l ) : 



E(r) = S e • exp{— b[(r/r e 



,1/n 



1]} 



where S(r) is the surface brightness at a radius r and S e is the characteristic value (i.e. the 
effective surface brightness), defined as the brightness at the effective radius r e . Parameter 
b is related to the Sersic index n and chosen equal to 1.9992n — 0.3271 so th at r e remains 



the projected radius enclosing half of the light in this component (ICiottilll99ll ) 
The disk component is an exponential profile of the form: 

S(r) = S • exp(-r/r d ) 

where Sq is the central surface brightness, and r d is the disk scale length. 



(2) 



Decomposition was performed based on this model, with a Gaussian PSF. A total of 
twelve parameters were adjusted in fitting the galaxy image and retrieved as output from our 
decomposition: the total luminosity L, the bulge fraction B/T, the bulge effective radius 
r e , the bulge ellipticity e, the bulge position angle (fib, the disk scale length r<f, the disk 
inclination angle i, the disk position angle (fid, the centroiding offsets dx and dy, the Sersic 
index n, and the residual sky background level db. The B/T which is defined as the fraction 
of the total flux in the bulge component has been extracted as a quantitative indicator of 
morphology. B/T = 1 corresponds to a pure bulge, while B/T = to a pure disk. With all 
the objects detected by SExtractor flagged except the galaxy of interest, the g— band mask 
images were adopted not only in g— band but also in r— band, throughout the decomposition 
fitting with GIM2D. Figure [3] shows example (science, mask, model, and residual) images of 
different morphological types in GIM2D fitting. The returned \ 2 values of g— and r— band 
fitting have the mean values of 1.14 and 1.18 with deviations of 0.26 and 0.33 respectively, 
representing rather convincing fitting results. 
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To check the results obtained, we compared the parameters B/T of the 154 galaxies 
estimated from g— and r— band in Figure HI The B/T values obtained from both bands 
agree well with small amount of deviation. Among the three most deviated sources, SDSS 
J161222. 61+525827.9 is an early-type galaxy with highly centralized surface brightness dis- 
tribution in bluer band. SDSS J151723. 30+593517.0 is a peculiar galaxy and UGC5888 is 
an irregular. We adopted B/T values derived from r— band images throughout the follow- 
ing investigation since the decomposition results in the two SDSS bands hold quite good 
agreement. 



3.2. Visual Classification 



The morphological types of revised Hubble sequence (jde Vaucouleurs fc de Vaucouleurs 



1964 ) of our sample galaxies were classified by visual inspection based on features like bulge 
ratios, the presence of spiral arms and/or bars, signs of interaction, multiple nuclei etc. 
All galaxies in our sample were classified into six morphological classes: T = 0(E or SO), 
l(Sa), 3(Sb), 5(Sc), 7(Sd), and 9(Irr). Notice that we assigned an additional class 10 cor- 
responding to galaxies with peculiar morphology possibly related to galactic interactions, 
mergers, etc. We performed the visual classification in both SDSS g— and r— band images. 
The classifications done independently by the four of us, agree in over 90% of the sample. 
The classifications were verified for 36 of those galaxies whose morphological types were 
found in the NASA/IPAC Extragalactic Database (NED) and were found to be accurate 
to AT = ±1. 

Considering the possible effects that inclinations and bar structures may cause, we 
further divided sample galaxies into 3 types: barred galaxies, edge-on galaxies, and the 
rest defined as general galaxies for further consideration. Note that when defining edge-on 
galaxies, we adopted a standard of GIM2D fitted disk inclination angle i > 70. Tabled] shows 
the numbers and corresponding fractions of different morphological types in our sample. 
In Table [2J we present the numbers of galaxies observed by different MIR bands in each 
morphological type. It can be drawn from Table 2 that the morphological fractions E/S0 
S(Sa-Sd): Irr of optical-8/im and optical-24/im samples are 0.47: 0.51: 0.014 and 0.44: 0.54 
0.016 respectively. Both are consistant with E(E/S0-S0): S(SOa- Sdm): Im of 0.40: 57 



0.014 obtained from an optically selected sample of SDSS galaxies (IFukugita et al.l 120071 ). 



http://nedwww.ipac.caltech.edu/index.html 
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3.3. Comparison 

A comparison between Hubble T type and the bulge-to-total ratio B/T has been carried 
out, as is shown in Figure [51 Except for the peculiar galaxies which exhibit a diversity of 
B/T, there is a smooth inclination in B/T along Hubble sequence from T = 9 to T = 0, i.e., 
for normal galaxies, except the peculiar, the late-types exhibit lower B/T than the early-type 
ones. Hence B/T does act as a reliable measure of morphological types for normal galaxies in 
our sample. In order to statistically examine the properties of our sample, galaxies have been 
divided into two morphological types: the early-type (E/SO) with T =0, and the late- type (Sa 



Sb, Sc , Sd, and Irr) with T from 1 to 9. Consequently, according to lSimien fc de Vaucouleurs 



(119861 ). such a classification roughly corresponds to the division of B/T = 0.4. 128 out of the 
144 normal sample objects coincide with their morphological types classified with T. The 
divisions are shown in Figure [5j Since as a whole, divisions with the two methods agree with 
each other, we have adopted T type to divide the sample galaxies into either the early-type 
or the late-type in the following statistics. 



4. RESULT AND DISCUSSION 

4.1. Color-Color Diagram 

The optical-MIR color-color diagrams are shown in the left panels of Figure EJ and there 
are anti-correlatio ns between optical and MIR colors, which is consistent with the result of 



Hogg et al.l (120051 ). Such a trend can be explained rather naturally: the optically blue color 
is related to the recent star formation, indicating the existence of notable amount of dust, 
and the MIR dust emissions can be shown in both 8fim and 24/xm bands. On the other 
hand, optically red galaxies are always too old to contain plenty of dust, and hence present 
weak MIR dust emissions. 

The right panels of Figure show the relationship between MIR colors. It can be 
detected that [3.6]-[4.5] colors in our sample are quite blue, with an average of around -0.52, 
indicating that there is no extremely active galactic nuclei (QSO etc.) in our sample, because 
both the emissions of 3.6fim and 4.5/im of normal galaxies are dominated by decreasing 
stellar continuum of the old stellar population. Yet these two bands of very active galaxies 
like quasars are dominated by the power law spectra of quasars and consequently should 
exhibit redder [3.6]-[4.5] colors. Both colors of [3.6]-[8] and [3. 6] -[24] can roughly characterize 
the relative strengths of the MIR dust emissions of PAHs and VSGs. The stronger the dust 
emissions, the redder the colors. Therefore, it is not surprised to find that most of the 
peculiar galaxies of our sample at the redder part of these panels because they contain large 
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amount of dust for their violent star formation (ISanders fc Mirabel! 119961 ) 



4.2. Estimation of the Stellar Mass 



T he luminosity of 3. Gu m band is often treated as a tracer of stellar component] Wu et al. 



2005a: Davoodi et al. 



20061. etc) as well as a test of validity of mass determination (IHancock et al 



20071 ). Furthermore IHancock et al.l (120071 ) has compared the mass of clumps in Arp 82 de- 
rived from R band fluxes and broadband colors against luminosities of 3.6/im, and has found 
strong correlation between them. Hereby the relationship between the stellar mass and the 
3.6/im luminosity of our sample galaxies is examined. Based on the optic al photome t ries a s 
described in § 12. 2\ we calculated the stellar mass of our sample, following iBell et al.l ( 120031 ): 



log( 



-0.4 x {M, r , AB - 4.67) + [a r + b r X (g - r) AB + 0.15] 



(3) 



where M^ab is the r— band absolute magnitude, (g — t)ab is the rest- frame color in th e 
A B magnitude sy stem. The coefficients a r and b r come from Table 7 of IBell et al.l (120031 ). 
A ISalpeterl (119551 ) stellar mass initial function (IMF) has been adopted with a = 2.35 and 
O.1M < M < 1OOM . The distribution of the stellar mass for our sample galaxies is 
presented in Figure [TJ in a range between 10 9 M Q and 1O 12 M , with the average mass around 
10 11 M o , as intermediate mass galaxies. 

A tight correlation of the stellar mass against the luminosity of 3.6/im is detected for our 



sampl e, as is shown in Figure [71 Our result is consistent with that found by IHancock et al. 



(120071 ) but with smaller scatters, thus it confirms the capability of the 3.6/im luminosity as a 
measure of the stella r mass of galaxies. Based on our 145 galaxies and 24 clumps in Arp 82 
(IHancock et al.ll2007l ). we fit the relation between the stellar mass and the 3.6/im luminosity, 
and obtain a nearly linear correlation: 



logi 



;i.34 ± 0.09) + (1.00 ± 0.01) x log 



uL u [3.6/im] . 



(4) 



As is pointed out by ICharlotl (119961 ) and iMadau et al.l (119981 ) , the mass-to-infrared light 
ratio is relatively insensitive to the star formation history, and remains very close to unity, 
independent of either galaxy colors or Hubble types. This relation is probably due to the fact 
that in the near-infrared, older stellar populations may dominate both the galaxy luminosities 
and the stellar masses. Therefore such a correlation provides a proxy way to estimate the 
stellar mass of galaxies directly from the integrated 3.6/im luminosity. 
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4.3. MIR Properties and Morphology 



Figure [8] shows the 8/xm and 24/zm dust luminosities as the function of different morpho- 
logical types in our sample. For normal galaxies, except the dwarfs, along either the Hubble 
T types or the B/T ratios, there exist obvious declinations of the MIR luminosities from the 
late-type to the early-type galaxies, especially showing a steep change around the division 
of the late-type and the early-type galaxies. All the peculiar galaxies show relatively higher 
MIR luminosities, independent of their B/T ratios. The phenomenon could be attributed 
to the fact that early-type galaxies are dominated by older population and are deficient of 
dust, resulting in lower MIR dust luminosities, while the late-type contain larger amount of 
young stars and more dust thus present stronger MIR emissions. Peculiar galaxies which 
are undergoing strong star forming ac tivities contain gr eat amount of dust and thereby show 
averagely high MIR luminosities. As IWu et al.l (j2005al ) has pointed out, both the 8/xm and 
24/im dust luminosities can be used as measures of SFRs of entire galaxies; therefore, the 
correlations between the MIR dust luminosities and morphological types reflect a consequent 
relationship be tween the galac ti c SFRs and morp hological types. This also confirms the pre- 
vious results of lSandagd (119861 ) . iKennicuttl (119981 ). etc. As for the six dwarf galaxies, because 
they all have low mass of around 1O 9 M (from the previo us mass determin ation), so they 
contain less dust and thus show lower MIR dust emissions (IHogg et al.ll2005l ). 



The MIR dust-to-star ratios are also plotted against different morphological types, in 
Figure [9j The prominent correlations are also detected between the MIR dust-to-star ratios 
and both the Hubble T types and the B/T ratios. Furthermore, such correlations seem 
to be more obvious than those between the MIR luminosities and morphological types. 
Since both the 8fim and 24/im dust luminosities possess correlations with SFRs for normal 
galaxies, and the 3.6/xm luminosity is a reliable tracer of stellar component, such ratios can 
be treated as the dust-to-star ratios. The correlations between the MIR dust-to-star ratios 
and the morpholo gical types could represent the distribution of SFRs per unit stellar mass 
(1 Wen et al.l 120071 ) over different morphologies. Contrary to the behavior in Figure [BJ the 
dwarf galaxies are roughly consistent with the other late-type galaxies within the error bars, 
but still present a little lower dust-to-star ratios. This could be explained with the fact that 
the gravitation potentials of these low mass galaxies may not be strong enough to retain as 
much dust and gas against the radiation fields as those of the normal-mass galaxies, and 
hereby all the dwarf galaxies show slightly lower MIR dust-to-star ratios. Almost all the 
peculiar galaxies present higher MIR dust-to-star ratios. 

The ratio of the 8/im dust luminosity to the 24/im dust luminosity is compared with 
different morphological types in Figure [101 Considering the 8/im dust emissions mainly 
represent emission of PAHs heated by B type stars while 24/im dust emissions stand for 
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emissions from hot dust mostly heated by O type stars (jPeeters et al.ll2004l ). the ratio can be 
treated as the estimation of ratio between these two components. In general, 8/im(dust)-to- 
24/zm(dust) ratios remain constant on average along the Hubble sequence, while the larger 
scatter in distribution of the early-type probably arises from existence of nuclear AGNs 
which destroy PAH emis sions presumably due to photodestruction of the PAH molecules 



by EUV/X-ray photons flGenzel et al.lll998 



Siebenmorgen et al.l I2004J ) or the outer diffuse 



PAH emissions heated by the older stars (jSauvage et al.l 120051 ). This can also be seen in 
later discussion in Figure [T21 where AGNs present lower 8/im(dust) luminosities. It should 



be noted that the two dwarf irregular galaxies exhibit lower 8Atm(dust)-to- 24jitm(dust) ratios 
than the most of the late-type do, possibly due to their low metalli cities ( Engelbracht et al. 
2005h . 



The above result indicates that for the normal galaxies except the dwarfs, the ratio of 
these two dust components does not vary much against morphology. Therefore the 8/zm(dust) 
luminosities are as capable to trace the galaxy SFRs as the 24/im(dust) luminosities as 
Wu et al.l (j2005al ) has pointed out. 



4.4. Statistics of MIR Properties 

In order to compare the statistical MIR properties of different types of sample galaxies, 
we performed all the following statistics on the sub-sample of 125 galaxies which have pho- 
tometric information for all the optical and MIR bands. Since from Figure [1] and K-S test 
results this sub-sample does not exhibit notable differences from the 154-galaxy sample, we 
suggest that the following statistical discussions are representative of all the sample galaxies. 

Distributions of the MIR luminosities and dust-to-star ratios for both the early-type and 
the late-type galaxies in the sub-sample are presented in Figure [TTl It is clear that, in general, 
the late-type galaxies exhibit quite different MIR properties to early-type ones, presenting 
a distinct and statistically higher MIR luminosities and dust-to-star ratios. The Gaussian 
fitting is performed on each distribution, and the lines crossing intersection points of the two 
Gaussian distributions are presented as the division of the two morphological types. In Part 
A of Table [3j the median values and the mean values together with scatters of distributions 
of the early-type and the late-type galaxies are listed. Probabilities that the distributions of 
these two types can match are all smaller than 8.8 x 10~ 5 , indicating that they present quite 
distinct properties. Therefore, both the MIR luminosities and dust-to-star ratios can be used 



2 with metallicity 12 + log(0/_ff) = 8.55 and 8.35 respectively, about 3/4 and 1/2 the solar value 
|Asplund et al.ll2004h . 
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to separate the early-type galaxies from the late-type ones. Table H] displays specific values of 
divisions and the reliability of such classifications. We define the reliability of classification 
as the fraction of galaxies from the subsample that are selected by Hubble T types. For 
example, out of the 117 normal galaxies which have images in all IRAC bands, the T type 
criterion selects 64 late-type galaxies and 53 early-types. 47 of the 64 late-type galaxies have 
logz/Ly(8/im)/L© ^ 8.91 and thus are consistently classified as late-type, while 47 of the 53 
early-types have logv 'L v (8 \xm) j 'L® ^ 8.91 and thus are classified as early-type. Therefore 
the logi/L„(8//m) / X = 8.91 maintains reliability of 73% for selecting late-type galaxies and 
89% for early-types. Correspondingly, log uL u (24:fim) / L Q = 8.30 gives a reliability of 70% for 
the late-type and 77% for the early-type; logz/L I/ (8/im)/^L I/ (3.6yum) = —1.15 can give a 88% 
for the late-type and 83% for the early-type; and loguL^Aixm) /uL v (3.6fj,m) = —1.45 can 
give a 84% for the late-type and 85% for the early-type; Hereby the MIR dust luminosities 
especially the MIR dust-to-star ratios can be effective tools to divide the early-type galaxies 
from the late-type ones. 

The comparisons between MIR properties of star-forming galaxies and AGNs (see §2.11) 
are shown in Figure [121 Generally, star-forming galaxies present stronger MIR emissions 
than those possess AGN activities. From Part B of Table [31 the star-forming galaxies and 
AGNs show the probability of matching each other in 8/im and 24fim dust luminosites 
and dust-to-star ratios with less than 2.1 x 10~ 3 , indicating quite different MIR properties, 
although not so distinct as these between the early-type and late-type ones. 

Some of the sample galaxies are edge-on galaxies and some are barred galaxies. Will 
the galaxies in edge-on view or with bars present different MIR properties? We compare the 
distributions of the MIR dust luminosities and the MIR dust-to-star ratios between edge-on 
galaxies, barred galaxies and normal galaxies in Figure [13] and Figure [HI The statistical 
results are listed in Part C and D of Table El with the probabilities in matching each by 
few to tenth percent, indicating rather resembling distributions. Thus, neither galaxies in 
edge-on view nor galaxies with bars present statistical differences from the other normal 
galaxies in both MIR properties. Furthermore, although based on limited data points, one 
can still find out in Figure [13] that the edge-on galaxies can also be classified into early-type 
and late-type ones with the division criteria described in Table [H Therefore, the inclination 
of galaxies and the existence of bars do not affect our previous results. 

We have also checked the distributions of total sample of 154 galaxies in the above 
parameters and type comparisons, and yielded quite similar results, therefore the selection 
of this sub-sample does not affect the reliability of our statistical results. 
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5. SUMMARY 

We investigated the correlations between the morphological types and the MIR prop- 
erties of a local optically flux-limited sample of galaxies selected from the spectroscopic 
catalogue of galaxies in SDSS-DR4, cross-correlated with the Lockman Hole, ELAIS-N1 and 
ELAIS-N2 of Spitzer SWIRE survey. Aperture photometry has been performed on all these 
galaxies in all optical and MIR bands. Morphological classifications have been performed 
by both visual inspection and the bulge-disk decomposition with GIM2D. Our major results 
are as follows: 

(1) The presented analysis clearly shows that the bulge-to-total ratio B/T obtained by the 
bulge-disk decomposition is proved a qualified quantitative measure of the Hubble T types. 
Galaxies with earlier morphological types possess larger bulge ratios while later type ones 
have more dominant disk structures. 

(2) The 3.6/xm luminosity presents a tight correlation with the stellar mass, and this provides 
us a new tool to estimate the stellar mass. The empirical formula to calculate the stellar 
mass by the 3.6/xm luminosity is given in Equation HI 

(3) Except for the dwarf galaxies and a few peculiar objects, the MIR dust luminosities of 
8fim and 24/im exhibit correlations with either Hubble T type or the bulge-to-total ratio 
B/T. Such correlations are much more obvious if we used the MIR dust-to-star ratios 
(either uL u [8fim(dust)]/i/L u [3.6^im] or vL u [24fim(dust)]/i>L u [3.6fim\) instead of the MIR 
luminosities. 

(4) The MIR dust luminosity ratios of uL 1/ [8fim(dust)]/vL u [2Afim(dust)] turn to be roughly 
constant against the morphological types, especially for the late-type galaxies. Therefore, 
on average, the 8/xm(dust) luminosity can as reliably measure the global SFRs of normal 
galaxies as the 24/xm(dust) luminosity can, regardless of the morphological types. 

(5) Distributions of the MIR dust luminosities and the MIR dust-to-star ratios of both 
the early-type and the late- type galaxies are very different. The late-type galaxies present 
higher MIR dust luminosities and MIR dust-to-star ratios than the early-type galaxies. The 
MIR dust luminosities of 8/xm and 24/im and especially the MIR dust-to-star ratios of 
uL u [8fMm(dust)]/ uL u [3.6^,m] and i>L v [24:fj,m(dust)]/ uL u [3.6fMm] can provide an effective tool 
to distinguish the late-type galaxies from the early-types. 

(6) The star-forming galaxies and AGNs also present different statistical MIR properties, 
with the former showing higher MIR luminosities and MIR dust-to-star ratios. The statistical 
results show that either galaxies in edge-on view or galaxies with bars do not present quite 
different MIR properties from other sample galaxies. 
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Fig. 1. — Distribution of sample galaxies' SDSS r— band Petrosian magnitudes, the redshifts, 
B-band absolute magnitudes and the stellar mass (deduced from § 14. 2p . Black histograms 
represent distribution of all the 154 sample galaxies, while the red for the 125 sample galaxies 
which have photometric information in all MIR bands. 
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Fig. 2. — BPT diagram plotted for 134 sources with emission-line detections. The solid curve 
illustrates the criteria based on the model of Kewley et al.(2001). Star-forming galaxies are 
labeled with magenta circles and AGNs with green asterisks. 
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Fig. 3. — Examples of the science, mask, model and residual SDSS r— band images for our sample, from top 
to bottom for galaxies corresponding to our morphological types: E/SO, Sa, Sb, Sc, Sd, and Irr. Elliptical 
isophotes for aperture photometry are overlaid in the first column. 
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Fig. 4. — The bulge-to-total ratios B/T in the SDSS r— and g— band are plotted against 
each other for the 154 galaxies modeled in both bands by fitting a Sersic profile to the bulge 
and an exponential to the disk. Error-bar representing 99% confidence are plotted. Note 
that different morphological types are plotted, with red filled circles for peculiar galaxies, 
green for general, blue and orchid representing edge-on and barred types respectively. The 
diagonal is plotted for comparison. 
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Fig. 5. — This plot exhibits the correlation between the bulge-to-total ratio B/T and T- 
type. Different colors of filled circles are defined as in Figure SJ The mean value of each 
type (including barred, edge-on, and general types) are overlaid in black diamonds, with 
standard deviation shown. Dwarf galaxies are overlaid with red diamonds. Divisions of 
normal galaxies with early- type and late-type are plotted, by T and B/T = 0.4, with peculiar 
galaxies excluded. The typical errors of sample galaxies are plotted in the upper right corner. 
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Fig. 6. — Color-color diagrams of our sample. The left panels show the derredened g-r color 
against MIR colors, while the right for MIR color-color diagram. Symbols are defined as in 
Figure [51 The typical errors are plotted in the lower right corners. 



-25 - 



12 r o Hancock, et al. 2007 



1 1 T 



10 r 

9 r 



O 



/ cro 



o 
o 



/ o 



/ o 



3. 

CD 



6/ 



4 r 

2 - 

7 
-2 L 



- "D "O 



O 



5 6 7 8 9 10 11 

Log 10 i/L„[3.6/xm]/L 



Fig. 7. — Upper panel: the 3.6/im luminosity vs. the stellar mass of our sample. Lower 
panel: the ratio between stellar mass and the 3.6//m luminosity vs. the 3.6/im lum i nosity . 
Symbols are the same defined as in Figure [5j Open circles are from lHancock et al.l (120071 ). 
corresponding to mass clumps in Arp 82. The typical random errors of log vL v [3.6/im] 
luminosity and logM* of these 169 data points are 0.09 and 0.01 respectively. 
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Fig. 8. — Distributions of photometric MIR luminosities for our sample, subdivided by 
different morphological types. Left panels: (upper) 8/xm dust luminosity and (lower) 24/xm 
dust luminosity vs. T; Right panels: (upper)8/im dust luminosity and (lower) 24/im dust 
luminosity vs. B/T. Symbols represent the same definitions as in Figure El The dot-dashed 
line represents the divisions of early-type and late-type. The typical errors are plotted in 
the lower right corners. 
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Fig. 9. — Ratio of 8/im and 24/xm dust-to-star ratios plotted against T and B/T. Symbols 
are the same as defined in Figure [U The typical errors are plotted in the lower right corners. 
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Fig. 10. — Distribution of ratios between 8/jm and 24/jm dust luminosities subdivided with 
different morphological types. Symbols are defined the same as in Figure [H The typical 
errors are plotted in the lower right corners. 
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Fig. 11. — Distributions of MIR luminosities and dust-to-star ratios for early-type and 
late-type galaxies of the 125 sub-sample. A: distribution of 8/xm(dust) luminosities; B: 
24/im(dust) luminosities; C: 8/im dust-to-star ratios; D: 24/im dust-to-star ratios. The 
black histogram represents distribution of early-types, and the red for late-types. The fitted 
Gaussian of each distribution is overlaid, with dividing lines through intersection points of 
Gaussian marked. 
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Fig. 12. — Distributions of MIR luminosities and dust-to-star ratios in Star-forming galaxies 
and galaxies with AGN activities (as in Figure fTTl) . The red histogram represents distribution 
of Star- forming galaxies, and the black for AGNs. Division lines are the same as in FigurefTTl 
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Fig. 13. — Distributions of MIR luminosities and dust-to-star ratios of edge-on and not edge- 
on galaxies (as in Figure [TTI) . The black histogram represents distribution of not edge-on 
galaxies, and the red for the edge-on ones. Division lines are the same as in Figure [TTJ 
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Fig. 14. — Distributions of MIR luminosities and dust-to-star ratios in barred and unbarred 
galaxies (as in Figure [TT]) . The red histogram represents distribution of barred spirals, and 
the black for the unbarred. Division lines are the same as in Figure [HJ 
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Table 1. The numbers of different morphological types of sample galaxies 





E/SO 


Sa 


Sb 


Sc 


Sd 


Irr 


Normal 


Pec 


General 


53 


13 


13 


11 


1 


2 


93 


7 


Barred 


6 


6 


8 


7 


1 





28 


2 


Edge-on 


7 


6 


5 


4 


1 





23 


1 


Total 


66 


25 


26 


22 


3 


2 


144 


10 


Fraction% 


42.9 


16.2 


16.9 


14.3 


1.9 


1.3 


93.5 


6.5 



Note. — Respective numbers and fractions of different morpholog- 
ical types and classes of 154 sample galaxies. Note that throughout 
the discussion in this paper, unless particularly stated, general, barred 
and edge-on types are referred to normal galaxies (from T=0 to T=9), 
while peculiar as an independent class. 
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Table 2. The numbers of different morphological types observed by each Spitzer MIR band 



Band 


E/SO 


Sa 


Sb 


Sc 


Sd 


Irr 


Pec 


Total 


3.6/mr 


61 


22 


26 


21 


3 


2 


10 


145 


4.5^m 


66 


24 


25 


21 


2 


2 


10 


150 


5.8/im 


61 


22 


26 


21 


3 


2 


10 


145 


8%m 


66 


24 


25 


21 


2 


2 


10 


150 


24/zm 


57 


22 


23 


22 


3 


2 


8 


137 


all bands 


53 


19 


20 


21 


2 


2 


8 


125 



Note. — The list shows respective numbers of different 
morphological types in each of the four IRAC bands and 
MIPS 24^m band, and also in the sub-sample whose 125 
galaxies have images in all MIR bands. 



Table 3. MIR Properties Distribution and Statistics 



Parameters 


Median 




Mean(scatter) 


Number 




D 


P 




(1) 


(2) 




(3) 




(4) 




(5) 


(6) 




< Part A > 


Early 


Late 


Early 


hate 


Early 


Late 








log vL v (&nm) /L Q 


8.30 


9.30 


o.zo^U.Uy J 


y.io^u.iuj 


53 


64 


0.72 


2.9 x 10" 


-6 


log vL u (24fim)/L 


7.76 


8.67 


( . / / (u.uy ) 


0.00(U.1U j 


53 


64 


0.63 


8.8 x 10" 


-5 


log vL v (8fim) 1 (3.6/zm) 


-1.55 


-0.22 


-1.4y(U.U ( J 


-U.oo(U.L) / j 


53 


64 


0.83 


4.7 x 10- 


-8 


log vL v (24^m) / ^L„(3.6^m) 


-2.06 


-0.81 


-z.UU(U.U8 j 


-U.yz(U.L) / j 


53 


64 


0.77 


6.9 x 10" 


-7 


< PartB > 


Star-forming 




otar-jorming 




Star-forming 


AGN 








log vL v (8jj,m)/LQ 


9.31 


8.78 


o nn/n i a \ 

o.yy(u.i4 J 


0. / Z(U.1U J 


51 


46 


0.38 


1.8 x 10" 


-3 


log vL v {2A^m)/ 'L Q 


8.75 


8.27 


o.4y^U.lo J 


O.IO^U.IUJ 


51 


46 


0.38 


2.1 x 10- 


-3 


log vL v (8/im) / vL^(3.6fj,m) 


-0.13 


-1.16 


-U.o ( ^U.IU j 


-i.u4^u.uy ) 


51 


46 


0.38 


1.8 x 10- 


-3 


log vL v (24^m) / vL v (S.&iim) 


-0.70 


-1.65 


-0.87(0.09) 


-1.60(0.10) 


51 


46 


0.38 


2.1 x 10" 


-3 


< PartC > 


Edge- on 


Not 


Edge-on 


Not 


Edge- on 


Not 








log vLv{&nm) 1 'Lq 


8.14 


8.81 


8.52(0.22) 


8.80(0.08) 


21 


96 


0.33 


0.05 




log vL v {2&nm) / Lq 


7.82 


8.31 


7.96(0.22) 


8.21(0.09) 


21 


96 


0.31 


0.08 




log vL v {%[im) /vL v (3.6 /xro) 


-1.26 


-0.75 


-0.94(0.19) 


-0.83(0.08) 


21 


96 


0.25 


0.25 




log vL v (2A:[ira) / Vi„(3.6^m) 


-1.83 


-1.37 


-1.51(0.18) 


-1.42(0.09) 


21 


96 


0.22 


0.40 




< PartD > 


Sarrerf [7 


nbarred 


Barred 


Unbarred 


Barred U 


nbarred 








log vL v {%nm) / L Q 


9.30 


9.44 


9.12(0.16) 


9.26(0.15) 


19 


31 


0.24 


0.51 




log isL v (24:fim)/LQ 


8.66 


8.78 


8.54(0.16) 


8.67(0.14) 


19 


31 


0.17 


0.86 




log i/L u (8jj l m)/i'L v (3.6fj,m) 


-0.15 


-0.22 


-0.26(0.12) 


-0.27(0.09) 


19 


31 


0.16 


0.91 




log uL v {2A^m) / vL v {3.%^Lm) 


-0.73 


-0.81 


-0.84(0.14) 


-0.86(0.08) 


19 


31 


0.21 


0.68 





Note. — Col. (2) are median values of the two compared types of sample galaxies. Col. (3) are mean values with scatters. Col. (4) 
are the respective numbers in statistics. Col. (5) are K-S test D (discrepancy) with Col. (6) the P (probability) that the two distributions 
match. The four parts correspond to: Part A are early-type and late-type sample galaxies; Part B are Star-forming and AGN samples; 
Part C are edge-on and not edge-on samples; Part D are are barred and unbarred spirals. 
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Table 4. Divisions of Early-type and Late-type Galaxies 





Parameters 


Divisions 


Reliability (L) 


Reliability(E) 




(1) 


(2) 


(3)' 


(4) 


log 


\vL„{%iim)/LQ 


8.91 


73% 


89% 


log 


;vL u (24[im)/L & 


8.30 


70% 


77% 


log 


; uL I/ {8jim)/vL^(3.6fim) 


-1.15 


88% 


83% 


log 


',vL v (2Aiim) 1 VL„(3.6^m) 


-1.45 


84% 


85% 



Note. — Col. (2) are the values corresponding to the division of the two 
types. Col. (3) are the reliability for selecting late- type galaxies, while Col. (4) 
are the reliability for selecting early- types. 



